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Abstract: We solve the Klein-Gordon equation for a scalar field, in the background ge¬ 
ometry of a dust cloud collapsing to form a black hole, everywhere in the (1+1) spacetime: 
that is, both inside and outside the event horizon and arbitrarily close to the curvature 
singularity. This allows us to determine the regularized stress tensor expectation value, 
everywhere in the appropriate quantum state (viz., the Unruh vacuum) of the field. We 
use this to study the behaviour of energy density and the flux measured in local inertial 
frames for the radially freely falling observer at any given event. Outside the black hole, 
energy density and flux lead to the standard results expected from the Hawking radiation 
emanating from the black hole, as the collapse proceeds. Inside the collapsing dust ball, 
the energy densities of both matter and scalar field diverge near the singularity in both 
(1+1) and (1+3) spacetime dimensions; but the energy density of the field dominates over 
that of classical matter. In the (1+3) dimensions, the total energy (of both scalar field and 
classical matter) inside a small spatial volume around the singularity is finite (and goes 
to zero as the size of the region goes to zero) but the total energy of the quantum field 
still dominates over that of the classical matter. Inside the event horizon, but outside the 
collapsing matter, freely falling observers find that the energy density and the flux diverge 
close to the singularity. In this region, even the integrated energy inside a small spatial 
volume enclosing the singularity diverges. This result holds in both (1+1) and (1+3) space- 
time dimensions with a milder divergence for the total energy inside a small region in (1+3) 
dimensions. These results suggest that the back-reaction effects are significant even in the 
region outside the matter but inside the event horizon, close to the singularity. 
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1 Introduction, Motivation and Summary of Results 


It is well known that, in the presence of a gravitationally collapsing structure forming a 
black hole and a quantum held in the Unruh vacuum state, an observer far away from the 
black hole will see a hux of thermal radiation at late times [1, 2], This result, which arises 
from the study of quantum helds in the curved spacetime, has led to several fascinating 
developments (see e.g. the textbooks and the reviews, [3-12]) in general relativity. While 
probing this result from different perspectives, quantum held theory in the region outside the 
black hole event horizon has been studied extensively in the literature. However, somewhat 
surprisingly, there has been much less emphasis in the study of quantum held theory inside 
the event horizon (for some earlier work, similar in spirit, see e.g., [13-18]). The purpose 
of this paper is to investigate quantum held theory inside the horizon in the context of a 
collapsing dust sphere in (1 + 1) spacetime. As we shall see, such a study leads to several 
curious and interesting results which should have their counterparts in (1+3) dimensions. 

In order to elaborate what is involved, let us consider the Penrose diagram in Fig. 1 
(top left) describing a gravitationally collapsing body. It is clear from the hgure that there 
are four distinct spacetime regions - marked A, B, C and D. Of these, region D — which 
is inside the collapsing body and outside the event horizon — is the least interesting one 
for our purposes. Even though the time dependence of the metric will lead to particle 
production in this region, we do not expect any universal behaviour here; the results will 
depend on the details of the collapse. Let us next consider region C which is outside both 
the collapsing body and the event horizon. This region is of primary importance — and 
has been extensively investigated in the literature — in connection with the black hole 
radiation. This is schematically illustrated in Fig. 1 (top right) by an outgoing null ray 
that straddles just outside the horizon and escapes to future null inhnity. The thermal 
nature of the black hole radiation arises essentially due to the exponential redshift suffered 
by this null ray as it travels from just outside the collapsing matter to future null inhnity. 
While this ray is inside the collapsing matter during part of its travel, the details of the 
collapse are sub-dominant to the effect of the exponential redshift at late times. We can 
investigate the black hole evaporation scenario vis-a-vis different kinds of observers in this 
region: like, e.g, asymptotic and non-asymptotic static observers, radial and inspiraling free- 
fallers, observers moving in circular orbits etc.; all these cases indeed have been studied in 
the literature (see, for some recent work. Refs. [19-21] which contain references to earlier 
papers). In this paper too, we will briehy discuss the physics in this region since recovering 
the standard results provides a ‘calibration test’ for our approach and calculations. 

But what we will concentrate on are the regions B and A which are inside the event 
horizon. (We have not found any extensive and systematic investigation of these regions in 
published literature which was one of the key motivations for this work.) The examination 
of the Penrose diagram in Fig. 1 reveals the following facts about these two regions. 

• Region B is inside the event horizon but outside the collapsing body. Being a vac¬ 
uum region in a spherically symmetric geometry, this region is indeed described by 
a Schwarzschild metric. But, if we use the standard (t, r, 6, (p) Schwarzschild coor¬ 
dinates, then r is like a time coordinate in this region due to the flip of signs in 
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the metric coefficients at r < 2M. Naively speaking, this makes the geometry “time 
dependent” (due to the dependence of geometry on r) in this region. Alternatively, 
one can describe this region using a coordinate system which is non-singular at the 
event horizon like, for e.g., the Kruskal coordinates, in which the line element takes 
the following form: 

ds^ = {-dT^ + dX^) + dLl ( 1 . 1 ) 

where, r is given as an implicit function of X and T via the relation: (r/2M — 

and the transverse line element is dL‘^ = r'^dO'^ + sin'^OdcjP'. 
Once again, the metric will be time dependent because of its dependence in the 
Kruskal time coordinate T. All these suggest that we expect a non-trivial dynamics 
for the quantum field in the region B. 

A typical null ray in this region, shown in Fig. 1 (bottom left), is trapped and hits the 
singularity in the region outside the collapsing matter in a finite time. (More precisely, 
it gets arbitrarily close to the singularity in the course of time.) It is instructive to 
compare the null ray in this region with the null ray in the region C mentioned 
above. (That is, we compare the figures at top right and bottom left of Fig. 1.) The 
two relevant null rays propagate straddling the event horizon just outside and just 
inside. While the outside ray reaches future null infinity (and plays a key role in the 
description of the black hole evaporation), it is not clear how the physical parameters 
vary along the ray which travels just inside the event horizon but gets arbitrarily close 
to the singularity at late times. This is important because the accumulation of energy 
density due to these modes of the vacuum which are trapped inside the event horizon 
can have important back reaction effects. One key aim of this work is to study the 
physical properties of the quantum field at the events along the null ray like the one 
shown in the bottom left figure of Fig. 1. 

• Region A is inside the collapsing body as well as inside the event horizon. Here too 
we can study the behaviour of the quantum field along the events in the path of the 
null ray shown in bottom right of the Fig. 1. This null ray also hits the singularity 
but inside the collapsing matter. The key difference between the events along the 
rays shown in bottom right and those in bottom left, is the following: In the former 
case, one could compare the energy density of the quantum field with that of the 
collapsing matter; We expect both to diverge as we approach the singularity inside 
the collapsing matter. The key question is to determine which one diverges faster in 
order to ascertain the effects of back reaction. But in the latter case (corresponding 
to the situation in bottom left figure), we have no matter energy density to compare 
with the energy density of quantum fields close to the singularity. Therefore, any 
divergence in the energy density of the quantum field will potentially have important 
consequences for the back reaction. 

We shall focus mostly on the regions A and B in this paper. We will use C mainly for 

“calibration"; that is, to arrange the quantum state in such a way that an asymptotic 
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r — 0 


r — 0 






Figure 1. The Penrose diagram in the top left figure illustrates the four spacetime regions A,B,C,D 
we are interested in. The other three figures describe the three null rays relevant to this work. See 
text for detailed discussion. 

observer sees the standard Hawking flux at late times. After choosing such a state, we will 
study the dynamical evolution of the quantum held and its regularized energy momentum 
tensor in the regions A and B. The region D is of comparatively less interest but we mention 
the results regarding this region in the appendix for the sake of completeness. 

Methodology of the paper. The Klein-Gordon equation for a massless scalar held can 
be solved exactly in (1 + 1) dimension, by exploiting the conformal invariance of the theory, 
which is a fairly standard procedure (see for a review [22]). The conformal invariance of 
the massless scalar held in (1 + 1) dimension (and the fact that any (1 + 1) dimensional 
metric can be written in conformally flat form) makes the relevant mode functions just 
plane waves. Therefore, the dynamics of a massless scalar held can be reduced essentially 
to ray tracing. As argued usually in the literature [7, 22], we expect the main conclusions 
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to carry forward to (1 + 3) dimensions, because the dominant contribution to the Hawking 
effect comes from s-waves even in (1 + 3) dimensions. 

Once the solution (j){x) to the Klein-Gordon equation is known (with suitable boundary 
conditions ensuring that one obtains the standard black hole evaporation in region C of the 
spacetime), our next task is to construct physically meaningful observables. There are 
two standard approaches which have been pursued in the literature in this context. One 
possibility is to introduce the particle detectors [23-26] in the spacetime moving on various 
trajectories. The particle content determined by the detector is then given essentially 
by the Fourier transform of the two-point correlation function of the held in the relevant 
quantum state. This leads to a Planckian spectrum for an asymptotic detector at rest in 
the Schwarzschild spacetime which agrees with the standard interpretation of black hole 
evaporation. Unfortunately, the response of the detector essentially measures the nature of 
vacuum huctuations and is sensitive to the history of the trajectory because it is dehned 
using an integral over the proper time. We cannot use the particle content determined 
by such a detector for estimating the back reaction effects of the quantum held on the 
spacetime. This is easily seen from the fact that a uniformly accelerated detector in hat 
spacetime will detect a thermal spectrum of particles but these “particles” do not back react 
on the spacetime in a generally covariant manner. 

Since our primary interest is to study the effect of quantum helds on the background 
geometry, we need to use a more covariant diagnostic. Such a diagnostic is provided by 
the (regularised) stress-energy tensor of the quantum held in the given quantum state. It 
is generally believed that, at least in the semi-classical regime, this regularized expectation 
value {Tab)reg Can be used as a source in Einstein’s equation to study the effects of back 
reaction. In particular, when {Tab)reg is comparable to the classical source of geometry Tab, 
we will expect back reaction effects to be signihcant. 

It should be noted that these two diagnostics for describing the quantum held — viz., 
the detector response or {Tab)reg — will, in general, give different results. Again, in the 
hat spacetime Minkowski vacuum, an accelerating detector will see a thermal spectrum but 
the regularized expectation value of stress-energy tensor will remain zero. For our purpose, 
it is clearly more meaningful to study (Tafe)reg rather than the detector response and we 
will concentrate on this study in this paper. However, in the last section of this paper, 
we will give the relevant results for the detector response, for the sake of comparison and 
completeness. 

The components of the stress-energy tensor depends on the choice of the coordinate 
system and hence could inherit the pathological characteristics of the coordinate system. 
It is therefore better to use physically well dehned scalar quantities that have an invariant 
meaning. The most natural scalar quantities at any event in spacetime can be constructed 
in the freely falling frame at that event which eliminates any acceleration effects. If we hll 
the spacetime by a suitable congruence of radially freely falling observers with four-velocity 
then we can construct [27-33] two useful scalar quantities, at any event in spacetime, 
by the following dehnition: 


U = {fab)u^u^-, T=-{fab)u^n^ 


( 1 . 2 ) 
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Figure 2. The collapse scenario in the Kruskal coordinates indicating the r = 0 and r = e surfaces. 
Of the two null rays which are marked, the ray 1 reaches arbitrarily close to the singularity in finite 
time while the ray 2 propagates to asymptotic infinity at late times. We will be interested in the 
energy density of the quantum field along the events in the path of ray 1, while ray 2 will be used 
for calibration. See text for detailed discussion. 

Here n“ is the normal in the radial direction (such that UaU'^ = 0), ^ is the energy density 
and T is the flux at that event as measured by a freely falling observer. 

As we said before, the most natural choice for is the four-velocity of the freely 
falling observers which are free of the acceleration effects. Further, we can All the entire 
spacetime with the freely falling observers, which allows a fairly uniform description of all 
the regions of the spacetime. But, in principle, one can also define the scalars in Eq. (1.2) 
corresponding to any observer with a given four-velocity and corresponding n“. This is 
relevant in regions C and D of the spacetime where one can also introduce static observers 
and compute the energy density and flux as measured by them. 

Before proceeding further one needs to clarify couple of issues related to this approach. 
First, note that the r = 0 singularity is a mathematically ill-defined event. Even though 
in the Kruskal coordinates in Fig. 2 (both in (1+1) dimension and in (1+3) dimension) it 
is drawn as a hyperbola, with distinct (T,X) coordinates along it, related by = 1 

one needs to be careful while dealing with the mathematically ill-defined nature of this 
event. (For example, consider the two events on which the two null rays in bottom-left 
and bottom-right figures of Fig. 1 hit the singularity. One could worry whether these two 
events are physically distinct.) This is particularly true in the (1+3) dimension in which 
the angular part vanishes at r = 0. To keep our discussion physically unambiguous, we will 
work throughout this paper with a surface infinitesimally close to r = 0 but non-singular. 
Since r 0 on this surface (see Fig. 2) any two separate points on it are indeed distinet. 
Thus by working on a r = e surface, with arbitrary small e these issues can be avoided. 

Second, the divergences in the densities in the (1 + 3) dimension can sometimes be 
spurious because the volume factor can become arbitrarily small. So one can have a situation 
in which U diverges but Uy/hd^x remains finite just because of geometrical considerations. 
So before we declare the energy density to be divergent in (1 + 3) dimension we need to 
ensure that the total energy contained inside a small volume itself diverges, rather than 
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just the energy density. This requires us to define a 3-volume element inside the event 
horizon. (One cannot, of course, use Airr'^dr because inside the event horizon, r = constant 
surfaces are spacelike.) When we use a freely-falling observer with the four-velocity to 
define the appropriate 3-volume measure to use is the one corresponding to the same 
observer in the synchronous coordinates. Using this measure and integrating U^/hd^x over 
a small volume around the singularity, we can determine the proper nature of divergence 
near r = 0. We will show later that the energy density U diverges in both regions A and B 
in (1+1) spacetime dimensions. However in (1+3) dimension the divergence in region A is 
compensated by the shrinking volume measure (thereby making the energy inside a small 
region around the singularity finite), while the divergence in region B persists, even though 
it is milder (being only compared to in the (1+1) dimensions). 

We will work with the Oppenheimer-Snyder model [34-40], which corresponds to col¬ 
lapsing dust (matter with zero pressure) that forms a black hole. In this case, the spacetime 
inside the collapsing matter is described by the closed Friedmann metric. To compute the 
(regularised) stress-energy tensor of the field in a particular quantum state, we use the tools 
of the standard conformal field theory techniques. (We expect the results to map to the 
s-wave sector of the scalar field in the collapsing background in (1 + 3) case which is the 
usual assumption made in literature; see e.g., ref. [7, 22]) This essentially implies that at 
any point of the spacetime, we have ingoing and outgoing (spherical) waves and hence the 
whole manifold can be coordinatized using these null rays. The incoming ray comes directly 
from J~ and the outgoing ray comes from J~ after a reflection from the vertical line in 
the Penrose diagrams, representing r = 0. As for the vacuum state, we shall work with 
the natural in-vacuum that is uniquely defined on J~. The relation of this state with the 
Unruh vacuum and its effects in the region C for the case of thin null shell collapse have 
been discussed in the refs. [19-21] and will not be repeated here. 

Summary of results. Our study leads to the following results, particularly in the region 
inside the the event horizon: 

• The energy density and flux, as measured by the radially in-falling observers, both on 
the inside and outside the dust sphere (but inside the horizon) diverge on approaching 
the singularity. (This feature was also observed previously for a null shell collapse 
[20]). Hence, the energy density for the scalar field can be arbitrarily large close to 
the singularity, even outside the dust sphere. The same results hold on events along 
the null rays, inside the event horizon (shown in the bottom figures of Fig. 1). For the 
ray reaching arbitrarily close to the singularity outside of the dust sphere (viz., the 
one shown in the bottom left of Fig. 1), both energy density and the flux diverge to 
arbitrarily high values as the ray approaches the singularity. Its implications remain 
to be ascertained but it could potentially affect the nature of the singularity (see e.g. 
[44, 45]). 

• Further, the total energy obtained by integrating the energy density over a small 

three volume near the singularity also diverges in both (1+1) dimension [as e“^] and 
in the (1+3) dimension [as This divergence might prohibit the formation of 
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singularity due to arbitrarily large back-reaction from the quantum field. It should 
be noted that whenever we probe the null rays or in-falling observers that hit the 
singularity, we can always work with spacetime events arbitrarily close to but distinct 
from the singularity, which have increasingly high curvature. 

• Let us consider the energy densities along the events in the rays shown in the top right 
and bottom left of Fig. 1. The spacetime events along these two rays, with one outside 
and one inside the event horizon, have approximately equal energy density when they 
are inside the matter. However, as we probe the events along the inside ray, which 
approaches singularity, the ratio of the energy densities on the outside events (which 
are approaching future null infinity) to the inside events (which are approaching the 
singularity) goes to zero due to arbitrarily high value of energy density, close to the 
singularity, on the latter null ray. 

• Similar behaviour is seen on the spacetime events along the null ray which are com¬ 
pletely inside the dust sphere. In this case as well, as we approach the singularity, 
we find a divergent energy density and flux. The ratio (U/p) of the energy density 
of the scalar field, to that of the dust also diverges as we approach the singularity in 
both the (1-1-1) and (1 -|- 3) spacetime dimensions. Hence the back-reaction due to 
the scalar field will affect the region inside the matter as well. 

• As regards the outside region, we find that the energy density and flux, measured by 
static observers and radially in-falling observers, exactly mimic the results obtained 
earlier in the literature in the context of a null shell collapse. This suggests that these 
results are generic and independent of the nature of collapse. The energy density 
at the spacetime events along the null ray moving forward to future null infinity is 
always finite and reaches the standard value (corresponding to Hawking evaporation) 
at late times. This is in agreement with the earlier investigations. 

The plan of the paper is as follows. We briefly review the matching of the interior and 
exterior parts of the spacetime for a collapsing dust ball in Sec. 2 and then introduce the 
double null coordinate spanning all of the spacetime. In Sec. 3, we consider the scalar field 
on this background and compute the components of the regularised stress-energy tensor 
for same. Subsequently, we compute the energy density and fluxes for various cases in 
Sec. 4. Introducing the radially infalling observers, we study the behaviour of the invariant 
observables U and T along the three null rays mentioned earlier in Sec. 5. Lastly, for the 
sake of completeness, we discuss the results regarding the detector response in terms of the 
effective temperature in Sec. 6. The last section contains the concluding remarks. 

2 The Gravitational Collapse Geometry 

In this section, we briefly review the junction conditions for matching the interior Friedmann 
universe to the exterior Schwarzschild spacetime.We will then introduce a useful global 
coordinate syatem (which we call the double null coordinates) and express both the interior 
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and exterior coordinates in terms of this double null coordinates in a conformally flat form, 
in the (1 + 1) sector of the spacetime. This allows us to determine the conformal factor of 
the metric from which we can calcualte the regularised stress-energy tensor. 


2.1 Junction Conditions 


We consider the collapse of a spherical region filled with pressure-free dust in the (H-3) 
spacetime dimensions. The inside region of the dust sphere is homogeneous and isotropic 
and the metric interior to the dust will be a closed (A; = 1) Friedmann model with the line 
element 

dsint = + a^iT)dx^ + a‘^{T)dL]_ ( 2 . 1 ) 

with r denoting proper time of the dust particles, comoving with = constant and = 
sin^ydfl^. A more convenient form for the line element can be obtained by transforming 
the proper time to conformal time r] via the relation: 


V = 



( 2 . 2 ) 


in terms of which the interior metric reduces to the following form: 


dsM = «^(d) (-V + dx^ + dLl) . 


(2.3) 


The Einstein equation for the interior Friedmann universe filled with dust can be solved in 
a parametric form leading to the results: 

a(d) = ^amax (1 + COST]) (2.4) 

Tiv) = ^dmax iv + sinry) (2.5) 

In these conformal coordinates the surface of the dust sphere is taken to be located at some 
value X = Xo and the collapse starts at jy = r = 0 and ends at ry = tt, r = {nl2)amax- The 
total energy contained within the dust sphere is constant and can be determined in terms 
of the quantity Umax as: 

o 3 

pa = constant = —Umax (2-6) 

OTT 

The exterior region, which is spherically symmetric and empty, is described by Schwarzschild 
metric but in a different set of coordinates. The spherical symmetry both inside and outside 
suggests that the angular coordinates for both the metric can be taken to be identical. Hence 
the outside line element has the following expression (we shall use the units with 2M = 1 
henceforth): 

= - ^1 - dt^ + (^1 - + dLl (2.7) 

where we have The above line element can also be written in terms of 

outgoing Eddington-Finkelstein coordinates v = t + r*^ with r* = r -|- In (r — 1), such that 
the line element reduces to 

dslxt ~ ~ ~ “^dvdr + dL\ (2.8) 
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In the outside coordinates the surface of the collapsing matter is characterized by r = r(r) 
and t = t{T), where, r and t are the Schwarzschild coordintes with r being the internal 
time coordinate. The above surface is also determined by x = Xo- Therefore the connecting 
equation for radial coordinate is [10, 38]: 


r (r/) = a{7]) sin xo = ^“max sin Xo (1 + cos rj) 


(2.9) 


Next we need to solve for the time coordinate, which is obtained by solving the following 
differential equation: 


1\ dt 

1-— = constant 

r J dr 


( 2 . 10 ) 


Matching the extrinsic curvature on both sides of the surface of the dust sphere we can fix 
the constant to be cos xo- This leads to the following differential equation for Schwarzschild 
time as: 


dt 2 d 

— O^max COS XO COS — T 
dr] 2 


®max COS XO COS 2 
Umax sin Xo COS^ 2-1 


( 2 . 11 ) 


Thus the time coordinate can be obtained from the above equation as: 


+ In 


( 2 . 12 ) 


The time corresponding to horizon crossing of the collapsing surface can be obtained by 
setting r = 1, which leads to, r]f{ = tt — 2xo- Note that as r/ —>■ rjH, t{rj) diverges; thus, 
for an outside observer, the surface of the imploding matter takes infinite time to reach the 
event horizon. The Umax and xo are not independent, since mass of the imploding dust ball 
is: 

dvr o dvr o q 1 q , , 

M = -^pr = -^PO' sm'^ Xo = -amax sm'=‘ xo (2.13) 

Thus with 2M = 1 we are left with the condition Qmax sin^ xo = 1- With these conditions we 
find the connection between the Eddington-Finkelstein coordinate v in the exterior region 
to that of interior region to be: 


v{ri) = jcotxo ( 1 + 
+ 2 In 


^max Sin 


r] + ^max sin Xo 


7 ] 7 ] 

cos - + cot Xo sm 


T] 

cos 2 


r] rj 

sm - + cot Xo cos - 


(2.14) 


Thus we have the matching condition of the outside Schwarzschild coordinates with the 
inside Friedmann coordinates. With this matching condition we now proceed to determine 
the double null coordinates for the entire spacetime region. 
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Figure 3. Penrose diagram showing the construction of double null coordinates. Each event x 
can be characterized by two null rays — an incoming ray from J~ which is labeled as and an 
outgoing ray, tracked backward to the vertical line r = 0 and then reflected towards J~ , giving the 
second null coordinate V~. 

2.2 Double Null Coordinates 

In this and subsequent sections we will first concentrate on the (1+1) spacetime dimensions 
by using dLj_ = 0 (which corresponds to d0 = 0 and dcp = 0) and quote the results for 1+3 
dimensions, whenever appropriate. For the purpose of calculating the vacuum expectation 
values for the regularised stress-energy tensor, it is useful to define and use a coordinate 
system made up of double null coordinates, constructed from the ingoing and outgoing null 
rays. Each event x (see Fig. 3) can be characterized by two null rays — an incoming ray 
from J~ ^ which is labeled as V~^ and an outgoing ray, tracked backward to vertical line 
r = 0 and then reflected to giving the second null coordinate V~. To obtain these 
null coordinates for the global spacetime we proceed as follows: In the interior Friedmann 
region we define two new coordinates which are both null 


(2.15) 


u = 7]-x; V = 7] +X 


Then the line element with the coordinates {V,U) in the interior Friedmann universe reduces 
to the following form: 


dsint = -a^{v)dUdV 


(2.16) 


In order to describe the exterior structure we define ingoing and outgoing Eddington- 
Finkelstein coordinates v and u both being null. The exterior line element reduces to: 



(2.17) 
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in the null coordinates {v,u). Then we can define the double null coordinates V~^ and V 
through the following relations related to the interior Friedmann universe as: 

V+ = A(V-xo) ( 2 . 18 ) 

V-=A{U-xo) (2.19) 


where we have introduced the function A[x) for future convenience: 

^max sin Xo 


A{x) = jcotxo (^1 + 
+ ^max sin Xo 


X + 2 In 


sin I + cot Xo cos | 


X X 

cos - +cotxosin-j 


X 

cos - 


( 2 . 20 ) 


which connect the interior coordinates (C/, V) to the global double null coordinates (1^^, V~). 

From the construction it is easy to verify the following results. The surface r = 0, which 
is in the FRW universe, has equation x = 0 or equivalently (from Eq. (2.15)) U = V. Thus 
from Eqs. (2.18) and (2.19) it is evident that = V~. Thus the two null coordinates 
coincide at r = 0. Also these null coordinates should be continuous on the surface. Then 
using the fact that V = ry + xo on the surface of the star from Eq. (2.18) we get, = v, 
which it should. For X 7^ 0; the value of V~^ and V~ are never equal for all ry. This implies 
that even though the singularity forms at r/ = vr, in the double null coordinates it is ‘spread 
out’. The value of V~ for the outer surface of the star to reach the singularity is obtained 
by using U = n — xo; leading to, V~ = A(tt — 2xo). By construction this must be equal 
to V~^ of the point, where the ray enters the dust sphere, which is A(r]). This helps us to 
identify, the time at which that particular null ray has entered the sphere as tt — 2xo. This, 
in turn, fixes the surface of the star completely. Note that this is just a corollary of a more 
general result: The null rays entering the sphere, getting reflected at r = 0 and exiting the 
sphere satisfy the relation, dry = dx- Thus the ray starts at x = XO) then goes to x = 0, 
again comes out at X = Xo. The difference in the coordinate values for rj between the point 
of entering and exit of the ray corresponds to Ary = 2xo (see Fig. 3). Hence in the above 
situation, the final value of ry is vr, and the entry value of rj has to be vr — 2 x 0 ) leading back 
to the previous result. 

Another set of relations connecting exterior coordinates (u.v) to the global null coor¬ 
dinates (F'*', V~) are given by: 


v = V+ 

u = B{U + xo) = A{U + xo) - amax sin xo (1 + cos{U + xo)) 

2 + Xo 


- 2 In 


^max Sin Xo COS 


- 1 


( 2 . 21 ) 

( 2 . 22 ) 


Thus the interior and exterior line elements can be written in the double null coordinate 
system (F"*", V~) in the following manner: 


dA 


= —a 


U + V 


2 / dA (M 

' dU dV 


dV+dV- 


ds 


2 _ 
ext 



dB/dU 

dA/dU 


dV+dV- 


(2.23) 

(2.24) 
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Hence this double null coordinate covers the full spacetime and brings the (1 + 1) sector of 
the spacetime to the conformally flat form. This is especially suited to evaluate the vacuum 
expectation value of the regularised stress-energy tensor, which is our next task. 


3 Regularised Stress-Energy Tensor 


We consider a minimally coupled, massless scalar field on the background geometry de¬ 
scribed by the line elements in Eqs. (2.23) and (2.24). In two dimension the dynamics 
of the geometry is encoded in the conformal factor which allows us to obtain the vacuum 
expectation value of the energy momentum tensor for the scalar field. For this we follow 
the standard procedure [22, 39] and use the following expressions given in terms of the 
conformal factor: 


{T++) 


1 

12^ 


Idle 

2~^ 


3 (d+cV 


(r__) 


1 


Idle 

2~C~ 


3 fd.ey 

4 v“^y 


{T+-) 


1 


d+d-e 

e 


d+ed-e 
~c ^ 


(3.1) 

(3.2) 

(3.3) 


In the above expressions we have introduced a short hand notation: the symbol ± stands for 
respectively. The detailed computation and the explicit expression for the regularised 
stress-energy tensor, using these relations is given in the Appendix A.l. 


4 Energy Density and Flux observed by Different Observers 

We shall now compute the energy density and flux at different events in the spacetime using 
the regularised stress tensor. As we described earlier our primary interest is in the regions 
B and A inside the event horizon where no static observers can exist. In these regions, we 
will study the energy density and the flux in the freely-falling frame. But before we do that, 
it is useful to consider the region C and see how the standard results are reproduced. In 
this region C (unlike in A and B) we can introduce both static and freely falling observers 
and study the flux and energy density as measured by both kinds of observers. 

4.1 Static Observers in region C 

We will start with the energy density and flux as measured by the static observers in region 
C. This will verify our procedure by leading to the standard result of the Hawking radiation 
with the Tolman redshifted temperature. 

A static observer stays outside the event horizon at some fixed radius r > 1. Since 
the observer is not following geodesic he/she has to accelerate (by firing rockets) so as to 
remain stationary at that radius. The velocity components for this static observer is given 


13 - 














by: 


(4.1) 


y+ = 


V- = 


r — 1 


r — 1 


cot xo — tan f \ cos^ f ^ 


V 2 


cot Xo + tan j ) 


(4.2) 


for a fixed r. The outward normal, determined from the condition naV°‘ = 0, has the 
following components: 


n+ = F+; n“ = -F” (4.3) 

Thus the energy density and flux can be computed using Eq. (1.2) for the static trajectories 
as a function of V~^. For static observers V~^ acts as the proper time along their trajectories 
and we plot both the energy density and flux as a function of V~^ in Fig. 4. From this figure 
we see that both energy density and flux shows similar behaviour at large radii. Both of 
them start from small positive values and finally saturate at the standard thermal spectrum 
values with the Tolman redshifted value for the temperature. 

The key results we want to verify for the static observers are, of course, the late time 
energy density and flux. The late time limit corresponds to C/ —)• vr — 3x0) under which 
both the energy density and the flux lead to the following expressions: 


U = 


7tT\ 


H 


T = 


12 

7rr2 


2 


r — 1 


12 r - 1 


(4.4a) 

(4.4b) 


We see that, in the asymptotic limit, r —)• 00 (corresponding to a static observer at a 
large distance) above expressions reduce to the standard Hawking energy density and flux. 
The redshift factor, as is well-known, diverges as the horizon r = 1 is approached. From 
Eq. (4.4a) it is evident that in the near horizon limit the energy density reaches a maximum 
and then decreases with decreasing r, eventually becoming negative. Thus, even in the 
Oppenheimer-Snyder dust model there is a region of negative energy density just outside 
the black hole horizon, just as in the case of a null shell collapse noticed earlier e.g., in Ref. 
[20]. However the flux is always positive and diverges at the horizon. 


4.2 Radially In-falling Observers 

Having reproduced the standard results in region C we now turn our attention to the 
phenomena inside the event horizon. Since there can be no static observers in this region, 
it is best to study these effects in the freely falling frame of the radially ingoing observer. 
To provide the complete picture and to maintain continuity, we will study the results for 
the freely falling observers in the entire spacetime. This can be done from two different 
perspectives: 

(1) We can compute the energy density and flux on the events along the trajectory of a 
given freely falling observer. This has a clear physical meaning of what such a freely 
falling observer will see, as she plunges into the singularity. 
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Figure 4. Variation of the energy density and the flux as observed by a static observer as a function 
of the proper time. The first figure shows the variation of the energy density with proper time . 
All the static observers at different radii will ultimately observe the standard result for the Hawking 
radiation with the temperature modified by the Tolman redshift factor. (Note that while plotting 
the graphs we have normalized them to their values at infinity, i.e. 7rT^/12 for convenience.) The 
same situation is depicted in the second figure for the flux. Both start with a small value and then 
rise rapidly, ultimately saturating at the appropriate values for the Hawking radiation with the 
redshifted temperature. 


(2) We can examine the energy density and flux on the events along specific null rays 
we described in Sec. 1 (see Fig. 1). This will tell us how the energy density and 
flux behaves on a selected set of events in the spacetime which are physically well- 
motivated. 

In the section, we shall present the results based on the approach (1), viz. along the 
trajectories of the observers. We shall describe the behaviour of the energy density and flux 
on the events along the specific null rays in Sec. 5. 

Since part of the spacetime is covered by the collapsing dust ball, there exist two types 
of radially in-falling observers. The first set of observers are (i) those who are comoving with 
the dust sphere and remain inside the collapsing matter, while the second set of observers are 
(ii) those who are always outside the dust sphere. The first set of observers who stay inside 
the dust sphere ultimately reach the singularity in a finite proper time. These observers 
are in a Friedmann universe. The observers in a region outside the collapsing matter also 
hit the singularity in finite time but they live in a Schwarzschild geometry. Thus these two 
observers are geometrically distinct even though both of them hit the singularity in a finite 
proper time. We will now consider these two sets of freely falling observers and use them 
to probe the regions D, A and B. 

4.2.1 Radially In-falling Observers: Inside regions D and A 

Let us start with the observers who are inside the dust sphere and are comoving with it. All 
these observers start at r/ = 0 when the dust sphere starts to collapse and reach ry = vr in 
a finite proper time. The trajectory of these observers are characterized by the conformal 
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time rj. For such an observer comoving with the dust sphere, the four velocity is given by: 


= 1 (dA dA\ 

^/2a(r?) \dV' dU) 

u _ 

“ ^\dA/dV'dA/dU) 

and the unit normal has the following expression: 


na 
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a 


a ( 1 

^Aw 

1 fdA 



dA\ 

W) 


(4.5) 

(4.6) 


(4.7) 

(4.8) 


The comoving observer is characterized by the value of y which remains fixed throughout 
the trajectory (taken as x) and hits the singularity as rj varies from 0 to vr. The variation 
of energy density for the comoving observer with the proper time r] along the trajectory is 
shown in Fig. 5. We note that the energy density for the outermost observer at x = Xo 
remains positive and diverges as r/ —>■ vr. While the energy density as measured by other 
observers remain finite while reaching a maximum near r] = tt and then becomes negative. 
Finally the energy densities for all the observers diverge in the ?? —)• vr limit. This result is 
shown in Fig. 5. 

The energy density diverges as the observer hits the singularity. The ratio of the energy 
density measured by the observer to that of the energy density of the collapsing matter also 
diverges. This divergence can be obtained from the leading order behaviour of U and p 
near r = 0, which can be represented, for (1+1) spacetime dimensions, as: (see Eq. (A.23) 
in the Appendix A.2.2) 




^ ^max ^ 

6vra^ ’ 


P 


( 1 ) 


^raax sin^ Xo ^ 
2a 


1 

7b ^ 


(4.9) 


where superscript (1) denote the values of respective quantities in (1+1) spacetime dimen¬ 
sions. Thus, close to the singularity, the energy density of the scalar field dominates over 
that of the dust sphere. 

We next consider the total energy within a small volume with linear dimension e in this 
(1+1) spacetime both for the scalar field (i?*-^^) and classical matter (Flb)j, This energy is 
given by: 


= juVh^dvA^ ~ = jp\/W)dvA) = e; ^ 1 (4.10) 

Thus even the total energy of the scalar field within a small volume diverges and is much 
larger compared to the total energy of the classical matter within that small volume. Energy 
in the quantum field dominates over the energy of the classical background. 


16 - 













Figure 5. The variation of the energy density as the collapse progresses for an observer inside the 
sphere and co-moving with the sphere. The energy density diverges as the conformal time rj tends 
to TT, the instant when the dust sphere collapses to a singularity. 


Generalization to (1+3) spacetime dimensions Let us now consider the generaliza¬ 
tion to (1+3) spacetime dimensions. Since we are considering s-wave approximation, the 
energy density for the scalar held in (1+1) spacetime dimensions can be related to that in 
(1+3) spacetime dimensions through the result [7, 22]: x (l/lvrr^) = where 

and are the energy densities in (1+1) and (1+3) spacetime dimensions respectively. 
So, inside the dust sphere, we have the following expressions for the energy densities of the 
scalar held and the dust: 


= 


^ ^max 

247r2a® sin^ Xo' 


,(3) 


3 0,max _ 

Svr ’ 


1 

/9(3) 


(4.11) 


where the superscript (3) denotes the values of the respective quantities in (1+3) spacetime 
dimensions. Hence the divergence inUjp still persists in (1+3) spacetime dimensions. 

As we have argued before, the really important measure is probably not the energy 
density but the total energy contained in a small volume of size e around the singularity 
in the e —>■ 0 limit. From Eq. (2.1) we arrive at the volume element to be: y/hd?x = 
sin^ X sill ^ dxdOdcj). The total energy inside a small volume can be found by integrating 
the energy density inside a sphere and is given by (with e —)• 0): 


T(3) = f U^^)x/md^x ^ e] = [ p^^'^VW)d^x ^ ^ = 4 (4-12) 

Jo Jo E'i e 

It vanishes for both components (as for the dust and as e for the held), where e represents 
radius of the small sphere around the singularity. However the total energy in the scalar 
held dominates over that in the classical background. This suggests that inside the dust 
sphere, the effect of back-reaction can not be neglected in (1+3) spacetime dimensions. (In 
normal units, the ratio of the energy densities in Eq. (4.12) will go as Lp/e^; so the effect 
is signihcant numerically at Planck scales, which is not unexpected.) 
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4.2.2 Radially In-falling Observers: Outside regions C and B 

The last set of observers we will study are the radially in-falling ones, but outside the dust 
sphere. They start from some initial radius and then follow a geodesic, ultimately hitting 
the singularity. They, however, live in the Schwarzschild region of spacetime throughout 
their life. Any such in-falling observer is parametrised by the energy associated with the 
geodesic, or — equivalently — by the initial radius from which her suicide mission starts. 
These two quantities are related via: 


E = 



(4.13) 


where rj stands for the initial radius from which the free-fall begins. From the geodesic 
equation we can determine the four velocity of these observers, which turns out to be: 


= 


r — 1 


E-\ - 


r — 1 


V- = 


r — 1 


E + i E^- 


r — 1 


cot xo — tan f \ cos^ f ^ 


\ 2 ) 


cot Xo + tan ^ ^ ^ j cos2 ^ ^ 


The components of the normal are determined by the condition UaV^ = 0 which 
the following choice for the outward normal: 


n+ = 1/+; 


n“ = -V- 


(4.14) 

(4.15) 


leads to 


(4.16) 


However we also require the evolution of V~^ as the observer proceeds towards the singu¬ 
larity, i.e. we need V~^ as a function of the observer’s conformal time r/. This can also be 
determined from the geodesic equation. For the radially in-falling trajectory the solution 
to the geodesic equation can be written as: 

r(r/) = ^ (1 + COST/) (4.17a) 

_ y ■ 

T = Ti + {t] + sliiti) (4.17b) 

In this case as well, the conformal time r] varies in the same range 0 < r] < tt. We also 
have the following initial conditions: r(rj) = rj. From these two relations we can determine 
dr/dri and, on using Eq. (4.14), we get the differential equation satisfied by V^{r]): 

dV+ ^[l-{l/ri)]‘''^cos(i|/2)-sm(i|/2) 

- l-{sscM.,/2)/r,} - 

The above equation can be integrated to obtain V~^ as a function of the time rj. Then 
both the energy density and flux can be obtained from this prescription. The behaviour of 
both the energy density and the flux are shown in Fig. 6. We note that both the energy 
density and the flux show identical behaviour. For rj < 0, when the sphere has not started 
collapsing, the energy density and flux vanish. (With our choice 2M = 1 the asymptotic 
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Figure 6. The first figure shows the energy density for the static observers at various radii at late 
times compared to the radially in-falling observer who crosses those static observers at different 
times. The second figure depicts the corresponding results for the flux. For static observers near 
the horizon the flux at late time diverges to infinity as the horizon is approached, while the energy 
density becomes negative. For radially in-falling observers nothing peculiar happens at r = 2M, 
while both U and F diverge as r —0. This feature was also observed earlier for the collapse of a 
null shell in [20]. 


value of Hawking flux would turn out to be: (l/1927r). To normalize the figures so that 
Hawking flux turns out be unity, the figures are drawn with a rescaling of the y-axis.) As 
the observer moves forward in proper time both the energy density and the flux arise and 
— as the event horizon is approached — both the energy density and the flux rise rapidly. 
However the values are finite at the event horizon located at r]H ~ 2.824. This suggests 
another interesting aspect to study, which is the comparison between the energy density 
and flux as measured by a radially in-falling observer and various static observers they 
encounter along the path. For the static observers the energy density and flux diverge as 
the event horizon is approached, while for the radially in-falling observer the energy density 
and flux remain finite at the horizon crossing. We give below the expressions for the energy 
density and flux as measured by these observers at the horizon: 

Uh = (^- - (4.19a) 

Th = vrT^ (4.19b) 

For a radially in-falling observer from infinity we have E = 1, for which the horizon crossing 
energy density and flux becomes Uh ~ 32Uoo and Eh ~ 24T'oo- (Note that these results 
were obtained earlier in the context of a null collapse in [20].) 

However the energy density and flux finally diverge as the singularity is approached. 
The behaviour of the energy density and its integral over a small volume has the following 
expressions near the singularity: (see Eq. (A.25) in Appendix A.2. 3) 

= J ~ ^ (4.20) 
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where superscript (1) denotes the expression of the quantities in the (1+1) spacetime dimen¬ 
sions. Thus both the energy density and the integrated energy diverges as the singularity 
is approached. Hence even in region outside the dust sphere the backreaction is important 
near the singularity. 


Generalization to (1+3) spacetime dimensions The energy density of the scalar 
held, as measured by the radially in-falling observers, in (1+3) spacetime dimensions can 
be obtained by dividing the (1+1) energy density by 47rr^. (This is a standard proce¬ 
dure adopted in the literature; see e.g., [7, 22]). This leads to the following approximate 
expression for the energy density near the singularity: 




7^2 1 

967r2 7-5 


(4.21) 


This result, at the face of it, shows that the back-reaction effects will be quite signihcant 
close to the singularity even outside the collapsing dust sphere. However, we need to 
ensure that the geometrical factor arising from the shrinking of the spatial volume does not 
over-compensate the divergence. To obtain the proper volume element appropriate for the 
radially in-falling observer, let us start with the metric in the synchronous coordinates for 
an observer in free-fall, from a large distance: 


ds'^ = —dr'^ H— dR^ + r^dVt^ 
r 


r = 




2/3 


(4.22) 


Thus for a T = constant surface the volume element turns out to be: 


VhdPx = sin 6 dRdOdcj) = -{R — t) sin OdRdOdcj) 


(4.23) 


Integrating this energy density over a sphere of small radius e we get the total energy to be 


re+T rj 2 f-e+T 


[i(fi-T)] 


10/3 


{R-'r)dR^ 


(4.24) 


which still diverges in the e —)• 0 limit but only as Thus, in the Schwarzschild 

spacetime region (outside the dust sphere), close to the singularity, the energy density due 
to scalar held tends to arbitrarily high value. The volume factor does help in (1 + 3) but 
not completely. 

Thus, on the whole, the results suggest that the back-reaction is important both in the 
outside Schwarzschild regime as well as inside the dust sphere. This has the potential of 
changing the geometrical structure near the singularity both inside and outside the matter, 
due to the backreaction. 

We will end this section with a few comments on how our results compare with those 
obtained in some previous attempts [46-49]. In most of these studies, approximate expres¬ 
sions for {Tab) have been obtained using a fourth order WKB expansion for the held modes 
to get unrenormalized (Tat) and then eliminating DeWitt-Schwinger counterterms [50] to 
get the renormalized value. All these approximate results (which includes both massless and 
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massive fields) for vacuum expectation value of the stress tensor has been obtained outside 
the event horizon. In [51] the approximate renormalized stress tensor in four-dimensional 
spacetime was obtained in the interior of the event horizon. Our results follow from the 
s-wave approximation [7] in which the dominant contribution to the Hawking effect comes 
from the monopole term (i = 0) in the multipole expansion. This is also well justified since 
we are using the radial observers to foliate our spacetime. 

As we are interested in the region near the singularity, we will consider dominant terms 
in the observables when the limit r —)• 0 is taken. If we calculate the energy density, i.e., 
{Tah)u°‘u^ for radially in-falling observer using the approximate stress energy tensor given 
in the ref. [51], it diverges as the singularity is approached (which has been pointed out in 
ref. [51] itself), but more importantly as ~ 1/r^. This is exactly the divergence we have 
obtained through our analysis as well. Hence, the energy density expressions obtained under 
s-wave approximation and energy density calculated using approximate renormalized stress 
energy tensor given in [51] have similar divergent behavior near the singularity. Moreover, 
the renormalized energy momentum tensor obtained in [51] includes curvature coupled 
scalar field as well. Hence the divergent nature of the energy density is a generic feature 
independent of coupling with curvature. 

Using perturbations around the Schwarzschild solution and treating the renormalized 
energy momentum tensor as a source for this perturbation, it is seen in ref. [51] that 
Kretschmann scalar diverges more rapidly for the perturbation compared to the classical 
Kretschmann scalar. This key result also follows and gets verified in our analysis. Along 
with these, a related fact that curvature for the perturbation grows more rapidly than the 
background geometry itself is also consistent with our results. 

5 Energy density and flux on events along speciflc null rays 

Combining all the results for the static and the radially in-falling observers we can describe 
the variation of the energy density and flux on events along the null rays introduced previ¬ 
ously. For the events along the null rays inside the horizon the best probes are the radially 
in-falling observers, as they pierce through every spacetime event in the inside region i.e. 
both in the regions A and B [see Fig. 1]. Thus each radially in-falling observer will intersect 
the null rays at one unique spacetime event. As the initial radius of the radially in-falling 
observer is varied, the trajectory will intersect the null ray at a different but still unique 
spacetime point. All these null rays in the outgoing mode have constant V~ and hence the 
only parameter that varies along the null rays is U+. For the events along the null rays 
outside the horizon we can determine the energy density and flux by using either the freely 
falling observers or the static observers. 

We summarize below all the results obtained for these three rays from our analysis. 

Events along the null ray just outside the horizon: Let us start the discussion by 
considering the behaviour of the energy density and flux on the events along the ray which 
straddles the horizon just outside and ultimately escapes to future null infinity 77'*'. (This 
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Figure 7. The variation of the energy density and flux on the events along a null ray as a function 
of time have been shown. Both the curves saturate at late time to the Hawking radiation values 
implying that all the events on do receive the standard Hawking flux. The rise of both the 
flux and energy density at lower values of V~^ corresponds to the near horizon behaviour of these 
invariant observables. See text for more discussion. 


is the ray shown in top-right figure of Fig. 1.) As we have discussed earlier, there is nothing 
peculiar happening in regions C and D, and hence we will not bother to discuss those 
regions. After its reflection at r = 0 the energy density rises, reaches a maxima and then 
drops back to the Hawking value as the asymptotic infinity is approached [see the left plot 
in Fig. 7]. Similar behaviour is exhibited by the flux, which also reaches a maxima and then 
decreases ultimately reaching the Hawking value as the asymptotic limit is approached. 

Hence we conclude that, on the events along the null ray just outside the horizon, 
nothing strange happens and the results reproduce the standard black hole radiation, known 
in the literature. 

Events along the null ray just inside the horizon: The second null ray passes through 
the events which are inside the event horizon (See the bottom-left figure in Fig. 1). Initially 
this ray shares the same geometry as the previous one and hence has almost the same energy 
density and flux. It also straddles the horizon but— being inside the horizon— its ultimate 
fate is sealed; it has to hit the singularity in finite time [see Fig. 1]. Thus even though 
the initial events along these two rays share similar physical conditions, the final regions 
of spacetime encountered by these two rays are very different. One ends at the future null 
infinity as described earlier, while the other one ends at the singularity. 

The energy density and flux in the present case can be calculated using the radially 
in-falling observers. For different radially in-falling observers (parametrized by the radius 
from which the free fall starts) the null ray cuts these observers at unique and distinct 
points. Thus calculating the value of the invariant observables at every point enables us 
to obtain their variation on the events along the null ray. The energy density along this 
null ray starts to rise and ultimately diverges as the singularity is approached [see the first 
figure of Fig. 8]. The same feature is also seen in the flux as measured along this null ray. 

Thus we can conclude that the final phase of the journey for this null ray shows sig- 
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Figure 8. The energy density and flux along null rays very near — but inside — the horizon. 
These rays come out of the dust sphere after the surface of the dust sphere has crossed the event 
horizon. Since all these null rays hit the singularity in a finite proper time the energy density and 
flux diverge. See text for detailed discussion. 


nificant differences compared to the one discussed earlier. For events along this ray close 
to the singularity, both the invariant observables diverge, in striking contrast to the events 
along the previous null ray, which ultimately leads to the Hawking flux [see Fig. 9]. Naively 
speaking this makes the back reaction effects quite signihcant close to the singularity. How¬ 
ever, as we said before, such a divergence in the energy density can be compensated by the 
shrinking 3-volume in a region close to the singularity. So we needed to compute the total 
energy inside a small volume before we can conclude about the divergence. This question 
has already been addressed in Sec. 4.2.2. The hnal conclusion to be drawn is that these 
quantities do exhibit a divergence even in (1+3) spacetime dimensions. 

Events along the null ray inside the dust sphere: This is the last situation we need 
to consider. This null ray stays mostly inside the dust sphere and also hits the singularity 
while it is still inside. On the events along this null ray we can make a direct comparison 
of the energy density of the quantum held with that of the dust sphere itself. As it moves 
towards the event horizon the energy density again rises, reaches a maxima and then diverges 
as the ray hits the singularity [see Fig. 1]. The same behaviour was also noted earlier for 
radially in-falling observers inside the dust sphere. As these observers are used to dehne 
the invariant observables along the null rays they exhibit similar behaviour. 

In this case as well the ratio Ujp diverges as the singularity is approached, i.e. the ratio 
tends to arbitrarily high value as the null ray reaches arbitrarily close to the singularity. The 
same result was also obtained earlier in the case of radially in-falling observers inside the dust 
sphere. Even though the above conclusions were drawn from (1+1) spacetime dimensions 
they carry forward to (1+3) spacetime dimensions as well, as shown in Sec. 4.2.1. Thus 
we conclude that the respective energy density for the scalar held and its ratio with dust 
energy density diverges in (1+1) as the singularity is approached. In (1+3), even though 
the energy for the scalar held does not diverge due to the volume factor in (1+3) spacetime 
dimensions, the ratio diverges. Thus the backreaction is important within the 
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Figure 9. The ratio of the two energy densities along the two null rays straddling the event horizon 
— one just one outside, while another just inside. As the inside ray approaches the singularity this 
ratio reduces to zero since energy density along the inside ray diverges. However for the rays near 
the event horizon — but outside — the energy density is finite, making the ratio vanish at late 
times. For two null rays with almost identical at the time of reflection at r = 0 the energy 
density turns out to have similar values; it can be seen from the figure that the ratio at the beginning 
is almost equal to unity. 


dust sphere as well. As we saw earlier, this is the case at events in region B close to the 
singularity, where the back reaction effects are significant. 


6 Effective Temperature measured by Detectors 


6.1 Effective Temperature as a measure of detector response 


We believe, for reasons described in Sec. 3, that the regularised stress-tensor expectation 
value is the appropriate quantity to study in our case. Neverthless, for the sake of com¬ 
pleteness and to compare our results with those obtained earlier in the literature, we will 
also consider the detector response in this Section. 

Let 7 (r) be the trajectory of an asymptotic stationary detector expressed in terms of its 
proper time r. Then, under certain conditions, one can associate [19, 41, 42] a temperature 
T_ with this detector by: 


r_ = 


1 


V- 

V- 


( 6 . 1 ) 
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Figure 10. Behaviour of U and at events along the null rays hitting the singularity while they 
are within the dust sphere. For these rays the value of the U coordinate is greater than the value of 
the U coordinate at which the dust ball hits the singularity. The null ray initially shows a growth 
as seen for the radially in-falling observers inside the dust sphere. After reaching a maximum the 
energy density ultimately decreases, finally it diverges as it hits the singularity. Identical divergence 
can also be seen for the flux. See text for detailed discussion. 


(One can also associate a second temperature T+ in an analogous fashion; but we will not 
need it for our discussion.) The usefulness of the above dehnition arises from the fact that it 
can be applied to non-stationary, non-asymptotic observers as well. It can be easily verihed 
[19] that T_ (r) indeed leads to the Tolman redshifted temperature [40] as measured by static 
Unruh-DeWitt detectors. The approximate constancy of T_ can be expressed through the 
adiabaticity condition [41, 42] which can be expressed as: 


V- = 


tI 


< 1 


( 6 . 2 ) 


This condition allows us to consider the response of the Unruh-DeWitt detectors in a 
straightforward way, by evaluating T_(r). Hence with every trajectory we can introduce 
a temperature T_(r) along with its adiabaticity parameter r/_. This setup will provide 
another handle on the quantum held theory of the scalar held in the collapsing background 
geometry. 


6.2 Static Detectors in region C 

For a static detector the effective temperature corresponding to V~ exists which, in the 
late time limit, ([/—)• vr — 3xo) reduces to: 

(the general expression is given in Eq. (A.27) in the appendix.) This is precisely the 
Tolman redshifted temperature, showing the validity and use of this effective temperature 
formalism. The variation of the effective temperature with V~^, proportional to the proper 
time of the static detector has been shown in Fig. 11. We have also illustrated the variation 
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Figure 11. The effective temperature T_ reaches the Tolman redshifted Hawking temperature Th 
at late F'*" (left figure). In the near horizon regime the effective temperature r_ is positive and 
diverges as the horizon is approached (right figure). 


of the effective temperature at late time with radii, showing the divergence as the horizon 
is approached. Incidentally, the detector temperature allows us to study another question 
which is of interest by itself. Consider a collapsing structure which — unlike the pressure- 
free dust studied so far — exhibits the following behaviour. The dust sphere starts to 
collapse from certain radius, continues to collapse untill it reaches r = 2M -|- e (with 
e <C 2M) and then stops collapsing due to, say, internal dynamics and becomes static. 
What kind of radiation will be detected by a static detector at large distances? Previous 
studies, based on quantum field theory [43], has shown that the collapsing body (i) will 
emit radiation closely resembling the Hawking effect during the collapsing phase and (ii) 
the effective temperature will drop down to zero at late times. It worthwhile to study this 
situation using the detector response. 

We will consider the detector response i.e. effective temperature measured by static 
detectors in this collapse scenario. The relevant collapsing geometry can be easily achieved 
by applying the same equations as given by Eq. (2.9) with rj being restricted to the range: 
0 < ?? < TT — 2xo ~ 2e. Here e is connected to the quantity Vstop — f's he. the radial 
distance from the horizon at which the collapse stops. This relation can be easily obtained 
as: Vstcrp — = 2 (cot xo) Hence, given the value of e, we can obtain the effective 

temperature at all times. We find it to be nonzero and approaching the saturated Tolman 
redshifted Hawking value during the collapse. While after the collapse stops there is no 
formation of trapped region and thus the effective temperature drops to zero. This situation 
is being depicted in Fig. 12 for two different choices of e. It turns out that as e decreases 
the effective temperature more closely resembles the Hawking temperature. As e —)• 0, we 
recover the usual Hawking evaporation result. 

Let us now return to our main theme. Having discussed the energy density, flux and 
effective temperature for static observers, we will now proceed to determine the same for 
radially in-falling observers, both inside the dust sphere and outside. 
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Figure 12. The variation of the effective temperature with F+ as the surface of the dust sphere 
approaches the horizon. If, for some reason, the collapse is halted at r = 2M + e (before forming 
the horizon)the effective temperature — which initially raises towards the standard Hawking value 
— drops to zero when the dust sphere becomes static. This situation is shown for two values of the 
final radii in the two plots. 
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Figure 13. The effective temperature T_ for radially in-falling detectors comoving with the dust 
sphere and remaining inside it. It is evident that the effective temperature diverges as the detector 
approaches the singularity. 


6.3 Radially In-falling Detectors: Inside regions D and A 

Let us first consider the effective temperature measured by the observers inside the dust 
sphere. As the spacetime can be globally mapped by the coordinates (y~^,V~), we can 
follow the same procedure as adopted before to calculate the effective temperature T_. For 
completeness we provide the expression for the effective temperature T_ along the trajectory 
of the radially in-falling detector in Eq. (A.28) in the Appendix. The effective temperature 
depends on the proper time r] along the trajectory of the radially in-falling observer and 
is plotted in Fig. 13. It becomes arbitrarily large as ry —?• vr i.e. as the radially in-falling 
detector approaches the singularity arbitrarily close. 
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Figure 14. The first figure shows the behaviour of T_, normalized to the Hawking value Th, as 
the radially in-falling observer approaches the singularity. The next figure describes the variation 
of the effective temperature with the energy gap of the detector in the asymptotic and near horizon 
regime. The near horizon value is greater than the asymptotic value for unbound trajectories i.e. 
trajectories with E >\. 


6.4 Radially In-falling Detectors: Outside regions C and B 

Let us now consider the case of Unruh-DeWitt detectors outside the dust sphere moving on 
a radially in-falling trajectory. The effective temperature they measure along their radially 
in-falling trajectory is specified by the energy E. This is a non-stationary phenomenon, 
since— as the detector approaches the singularity the local curvature grows rapidly. Thus 
this trajectory allows us to probe the time dependent Hawking temperature in detail. The 
expression is given in Eq. (A.29) in the Appendix. 

Here we are interested in two limits: the asymptotic one and the near horizon limit. In 
the asymptotic limit we have il = 0 and (r — l)/r = 1, such that the effective temperature 
reduces to: 

jnAsym ^ (6.4) 

which is consistent with what we expect. For a radially in-falling observer who starts 
her journey from spatial infinity has energy E = 1. Then this observer will detect a 
temperature = Th- While for observers with E ^ 1, the asymptotic temperature 

would be different from the Hawking value, as obtained earlier in the context of a null 
collapse in [19]. In this regime the adiabatic parameter is negligibly small. 

On the other hand, at horizon crossing, we have the following expressions: r = 
—E, = (1/2E) and r = —(1/2). Then we obtain the following expression: 

V~ 2u — 4f^(r +1) — 4r(—f + rv) 

^— = - = iim-- - - — = —2E 6.5 

y- 2u 4rr 

Thus the near horizon effective temperature is given by : 

Thus the effective temperature is not only non zero, but for states with E >l,\t exceeds the 
Hawking value. For an Unruh-DeWitt detector dropped from infinity we have E = 1 and 
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thus the detector will perceive four times the Hawking temperature at the horizon crossing. 
(This result was obtained earlier in Ref. [19, 41]). However at this stage the adiabaticity 
parameter r/_ has large value, and hence the interpretation of T_ as the temperature has 
some ambiguity, and one needs to consider progressively larger frequency modes. 

The behaviour of T_ has been plotted against the conformal time rj for the detector in 
Fig. 14. The effective temperature remains finite at the horizon crossing rjn ~ 2.82, while 
diverges later as the detector hits the singularity. We have also plotted the behavior of 
and as a function of the energy of the in-falling detector. We note that for 

unbound trajectories {E > 1), the horizon crossing temperature is always greater than the 
asymptotic value. 

7 Conclusion 

While the quantum field theory outside the black hole event horizon is well studied in the 
literature, the corresponding issues in the region inside the event horizon have not attracted 
sufficient attention. This was the key motivation for this paper. We have considered a 
massless scalar field in the background geometry of a collapsing dust ball (Oppenheimer- 
Snyder model) in the in-vacuum state defined on . We focussed on the (regularised) 
stress-energy tensor in this vacuum state as a physically relevant diagnostic of the quantum 
effects. The stress-energy tensor, in turn, leads to two scalar observables - energy density 
and flux in the normal direction - related to an observer at any point in the spacetime. 
The key results obtained through this analysis are summarized below: 

• Let us first start with the radially in-falling observers inside the dust sphere and the 
invariant observables they measure along their trajectory. 

1. In (1+1) spacetime dimensions the radially in-falling observers inside the dust 
sphere find that the corresponding energy density and ffux becomes arbitrar¬ 
ily large as they approach the singularity (arbitrarily close) in a finite proper 
time. Even though the energy density of the dust ball (p) itself diverges as the 
singularity is approached, scalar field energy density {U) diverges faster. This 
can be seen from the ratio lAj which diverges as e“^. The same conclusion 
holds for the total energy within an infinitesimal volume as well. Hence in (1+1) 
spacetime dimensions the scalar field dominates over the classical source. 

2. In (1+3) spacetime dimensions also the energy densities of both the scalar field 
and classical matter becomes arbitrarily large close to the singularity. Just as in 
(1+1), in (1+3) as well the scalar field energy density dominates over the classical 
source. When integrated over a small volume near the singularity energy in the 
scalar field goes as e, while that in the classical matter goes as e^. So, even 
though both of them vanish, in the e —)• 0 limit, the classical source vanishes 
more rapidly. Thus even in terms of the energy the scalar field dominates over 
the energy in classical background. Thus backreaction is important inside the 
dust sphere. 
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Figure 15. Penrose diagrams (along with the associated graphs) illustrate how the flux along 
the null rays vary. In the top-left figure we consider a ray straddling the event horizon but remains 
outside. The flux in the near horizon regime first shows an increment, but soon tends to the 
Hawking value as the ray approaches the future null infinity. The behaviour of flux along the null 
ray straddling the event horizon just inside, but remaining outside matter at late times is shown in 
the top-right figure. This null ray hits the singularity in a finite proper time and the flux diverges. 
The bottom figure shows the divergent nature of the flux for null rays inside the dust sphere. The 
flux diverges near the singularity as expected. 


• The second class of observers we consider are the radially in-falling observers outside 
the dust sphere. For them, we arrive at the following conclusions for the invariant 
observables: 
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1. In (1+1) spacetime dimensions, radially in-falling observers outside the dust 
sphere hnd that the energy density of the scalar held becomes arbitrarily large 
near the singularity. The same is also true for huxes measured by these radially 
in-falling observers near the singularity. This conclusion holds even after inte¬ 
gration over a small volume around the singularity. Since there is no classical 
matter present in this region this divergent energy of the scalar held has the 
potential to alter the singularity structure. 

2. The corresponding situation in the (1+3) spacetime dimensions is similar and 
both the energy density and hux diverge as the singularity is approached. Even 
after incorporating the effect of shrinking volume element d^x\/h, the divergence 
persists, but becomes milder (e“^/^), compared to the divergence (e“^) in (1+1) 
spacetime dimensions. The divergence in the energy density of the scalar held 
can prohibit the formation of singularity itself due to existence of back-reaction 
near the singularity. 

Thus, the energy density in the scalar held dominates over the classical background 
energy density making the backreaction important both inside and outside the dust 
sphere. All these conclusions hold for energy density and hux as measured along the 
events in the null rays hitting the singularity either inside the dust sphere or outside. 

• The radially in-falling observers outside the dust sphere observe that there is nothing 
peculiar at the event horizon; neither the energy density nor the hux diverges there. 

• The study of the events along the null ray straddling on the outside of the event 
horizon leads to the standard Hawking energy density and hux at late times, while 
the events along the ray inside the horizon sees an increase of both energy density 
and hux to arbitrarily high values as it approaches near the singularity. Even though 
derived in (1+1) spacetime dimensions, the same results hold in (1+3) spacetime 
dimensions as well, with milder divergences. 

In addition to the above main results, we have also obtained several results for the 
observers on the outside of the event horizon which agree with the previous studies and 
expectations. Also, for completeness, we have studied the effective temperature formalism 
as well and observed similar effects of divergence at the singularity in this case as well. 
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A Appendix: Details of the calculations 


In this Appendix, we present the steps for deriving various results mentioned in the text. 


A.l Stress-Energy Tensor: Explicit Derivation 
A. 1.1 Exterior Region 


The various derivatives of the conformal factor in the exterior region have the following 
expressions: 


r - 1 dB/dU_ p 2 ^ _ (r- 1)(3 - 2r) dB/dU 
2r3 dA/dU' + “ ^ dA/dU 


(A.l) 


d-C 


r-1 /dB/dU\ 
r ~ \dA/dU) 


r-1 (dB/dUV 
2r3 \dA/dU) 


(A.2) 


die 


3r-ldB/dU fdB/dU\ r-l^^ fdB/dU\ (3-2r)(r-l) 
2 r3 dA/dU^~\dA/du)^ r “ \dA/du)^ 

„ r-1. (dB/dU\ (dB/dUV {3-2r){r-l) 


fdB/duV 

\dA/dU ) 
(A.3) 

(A.l) 


With the following expressions for energy momentum tensor: 
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these relations can be simplihed to arrive at, 
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(r__) = — 
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In arriving at the above relations we have used the following expressions for the various 
derivatives dA/dU and dB/dU are respectively: 
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A.1.2 Interior Region 

For the interior region the various derivatives of the conformal factors are: 
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Then the components of the stress energy tensor in the inside region are: 
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where various derivatives of the quantity A are given in Sec. A.1.1. 


A.2 Energy Density and Flux For Various Observers 
A.2.1 Static Observer 

Below we provide the full expressions for energy density and flux calculated for static 
observer: 
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and the expression for flux turns out to be: 
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A.2.2 Radially In-falling Observers: Inside 

The energy density for radially in-falling observer has the following expression: 
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A.2.3 Radially In-falling Observers: Outside 

For radially in-falling observer outside the dust ball has the following expression for energy 
density: 
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and the flux has the following expression: 
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A.3 Effective temperature for various observers 

In this section we present the full expression for effective temperature measured by detectors 
in various trajectories. For static observer the effective temperature turns out to be: 
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For radially in-falling observes inside the dust sphere the effective temperature turns out 
to be: 
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Finally, for radially in-falling observers in the Schwarzschild spacetime effective temperature 
takes the following expression: 
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